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Abstract 

The local persistent current in two dimensional strongly interacting systems is 

investigated. As the interaction strength is enhanced the current in the sample 

undergoes a transition from diffusive to ordered flow. The strong interacting 

flow has the characteristics of a plastic flow through dislocations in the pinned 

charge density wave which develops in the system at low densities. 
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The behavior of interacting electrons in random potentials has drawn much interest 
because of its relevance to many different phenomena such as persistent currents the 
two-dimensional MIT (metal-insulator transition) ||, and the charging spectra of quantum 
dots It is generally believed that the large amplitudes of persistent currents measured 

for mesoscopic rings in the diffusive regime |5],|7| are the result of the suppression of the 
influence of disorder by interactions, although no quantitative calculation for a realistic 
model of a 3D metallic system has been forthcoming. 

Much recent work has concentrated on the influence of interactions on 2DEG (2D Elec- 
tron Gas). It is speculated that the experimentally observed 2D MIT is due to e-e 
(electron-electron) interactions ||, while indications of the influence of e-e interactions in 
the addition spectrum of disordered quantum dots are mounting H[|. For typical 2DEG 
devices the density is rather low (n < 3 x 10 11 cm~ 1 ) which corresponds to a ratio between 
the the Coulomb and Fermi energy r s = e 2 /vf > 1 (where Vf is the Fermi velocity) 0. In 
this region the simple RPA (random phase approximation) no longer holds, and correlations 
play an important role ||. 

Persistent currents for a clean 2DEG ring has been recently measured and shown to be 
of the expected magnitude fill. No measurements have yet been performed in the diffusive 



regime, but numerical calculations indicate that the persistent current will be significantly 
enhanced (compared to the non-interacting value) by the e-e interactions even for spinless 
electrons [JTT| . 



In this paper we would like to clarify the mechanism of this enhancement in 2DEG rings. 
A useful clue may be obtained from the study of the distribution of the current in different 
realizations of disorder [12] . While in the non-interacting case the distribution is broad with 



almost the same probability for dia or paramagnetic currents leading to small paramagnetic 
average currents, in the interacting case the distribution is predominately paramagnetic 
leading to large paramagnetic average currents. In this paper we shall show that the local 
persistent current undergoes a transition from a diffusive behavior for the non-interacting 
regime to a strongly correlated plastic flow through a small number of channels for stronger 



e-e interactions. This leads to a change in the local distribution of the current and to the 
enhancement of the total average current. 

Thus, there are two distinct mechanisms for the enhancement of the persistent current 
in 2DEG systems. One which is relevant to the high density weak interaction regime which 
is described analytically by the diagrammatic perturbation treatment [I3"|,I4"|, or by the 



Hartree-Fock picture fi~5"l , and treated numerically by various Hartree-Fock approximations 



16| -|T9|. The other is the low density strongly interacting regime which is characterized by 



the appearance of correlations in the electronic density and plastic flow of the persistent 
current. 

The local persistent current of an interacting 2D cylinder of circumference L x and height 
L y threaded by a flux $ was calculated. This system is known to show a large enhancement 



of the persistent current in the diffusive regime [11]. The Hamiltonian is given by: 



U t tx m i tt a k,j a k,j a i tP a i,p , 1N 

H = 2^£k,jalja k j -i^ J kJ + J% d + U ^ -±- — j£— (1) 



k,j k,j k,j>l,p 

where 



rkj - n,. 



J k,j = ~iV exp(i$s/ L x )al j+1 a kJ - h.c, 

J h = - iVa k+i,j a k,j - h.c, (2) 

and a k j is the fermionic creation operator, e& )3 - is the energy of a site (k,j), which is chosen 
randomly between — W/2 and W/2 with uniform probability, V is a constant hopping matrix 
element and s is the lattice constant. The distance \ fhj — fi tP \ = (min{(/c — I) 2 , (L X J s — (k — 
I)) 2 } + min{(j — p) 2 , (L y /s — (j — p)) 2 }) 1 ^ 2 . The interaction term represents a Coulomb 
interaction between electrons confined to a 2D cylinder embedded in a 3D space with U = 
e 2 /s, 

We consider a 4 x 4 lattice with m — 16 sites and n = 8 electrons. The many-particle 
Hamiltonian then may be represented by a 12870 x 12870 matrix, which is exactly diago- 
nalized. The many-particle ground state |^($)) for $ = 7r/2 is calculated for 500 different 
realizations of disorder in the diffusive regime (W = 8V WW) for several values of interaction 
U. The local persistent current 



it d m = {*m\Jkj\*{*)), ( 3 ) 

(where a — x,y) is calculated for each realization. 

Local currents in typical realizations are plotted in Fig. [I]. It can be seen that the 
current for the non-interacting case U = is diffusive, the current flows in all directions, 
there are current loops and there is no obvious long range correlations. On the other hand 



for U = 10V, which corresponds to r s ~ ^Jit/2(U/AV) ~ 3 a value which could be easily 
obtained in contemporary 2DEG devices, the behavior is totally different. The persistent 
current flows through a few channels in the sample, there are no currents flowing in the 
opposite directions, no closed loops, and obvious correlations. This situation reminds us of 
plastic flow in lattices, for example the flow of magnetic vertexes in driven Abrikosov lattice 
pinned by random pinning centers JHJ . 



Indeed for interaction strength corresponding to r s > 1 the system described by the 
Hamiltonian in Eq. ([!]) is known to exhibit short range density correlations which develop 
into a CDW (charge density wave) for stronger interactions ||. The charge density and the 
persistent current for some typical realization are shown in Fig. ^. Most of the current flows 
in channels which correspond to dislocations in the pinned charge density wave clearly seen 
in the density plot of the different realizations. This is a feature which is characteristic of 
plastic flows. 

A more quantitative measure for the change that the persistent current in the sample 
undergoes due to the e-e interactions is given by the correlation between the currents at 
different locations. This correlation may be formulated in the following way: 

1 ' ; ((^($)) 2 > - (i x {®)) 2 

where 

((i a m N ) = m d m N ) ( 5 ) 

and (. . .) denotes average over different realizations of disorder. The average current (J a ($)) 



as well as the typical current J ((/ a ($)) 2 ) on a bond are plotted in the inset of Fig. ^|, while 



the current per bond distribution P(I X ) for different values of interaction strength are plotted 



in the x direction increases up to U = 10V and then decreases, while the typical current is 
somewhat enhanced up to U = 10V and then decreases. Of course the average current in 
the y direction is zero and the typical current is somewhat enhanced in the weak interaction 



regime (r s < 1 corresponding to U < y32/7r~3, but suppressed for higher values of 
interaction. 

The average current is strongly enhanced in the regime of r s < 3. It is clear from the 
typical bond current dependence on U and from the distribution that the enhancement is 
not due to an enhancement of the typical current in a certain realization but due to the fact 
that in the interacting case the local current in the x direction for almost all realizations is 
paramagnetic while in the non-interacting case there is an almost equal probability of the 
current being para or diamagnetic. 

The correlation between the persistent currents at neighboring bonds as function of the 
interaction strength is seen in the plot of C(r = 1, $ = 7r/2) presented in Fig. While in the 
weak interaction RPA regime there is no enhancement of the local current correlations due to 
the increase in the interaction strength, there is a strong enhancement of these correlations 
in the strong interaction low density regime r s > 1. Thus, the impression one gets from 
the current plots shown in Figs. [I] and of a transition in the current characteristics of the 
sample, is confirmed. 

The reason for the enhancement of the average total current becomes now quite clear. 
While in the weak interaction RPA regime there is some enhancement in the average current 
due to interactions, this enhancement is not related to correlations in the local current. On 
the other hand, in the strongly interacting low density regime the enhancement is connected 
to an ordering of the local persistent current. The local persistent current in the strongly in- 
teracting case must flow along dislocation lines. There is only a small probability of forming 
closed loops and changing directions resulting in the fact that the current in most realiza- 
tions flows in the same direction and the distribution of the current is almost exclusively 



in Fig. [| As in the case of the total average current shown in Ref. |TTJ the average current 
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paramagnetic. 

In conclusion, strong interactions (or low densities) impose a significant modification in 
the nature of the local persistent current in a 2D sample. The current is transformed from 
a diffusive current to a plastic flow along dislocations of the emerging pinned CDW of the 
system. This transformation causes a strong enhancement in the total persistent current 
of the system since there is no suppression of the average current due to the fact that in 
different realizations current might flow in arbitrary directions. 

We are grateful to the Israeli Academy of Sciences and Humanities research center 
"Strongly Interacting Electrons in Restricted Geometries " for financial support. 
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FIG. 3. The distribution of the bond current in the x direction for different values of U. Inset: 
average (full line) and typical (dashed line) bond currents in the x direction and typical (dotted 
line) bond currents in the y direction in units of the hopping matrix V. 
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FIG. 4. The bond current correlation in the x direction as function of the interaction strength 



U. 
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